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ABSTRACT 



Context. The star y Cep is known as a single-lined spectroscopic triple system at a distance of 13.8 pc, composed of a Kl III-IV primary star 
with V = 3.2 mag, a stellar-mass companion in a 66-67 year orbit (Torres 2006), and a substellar companion with M,, sin; = 1.7 Mji,p that is 
I ' most likely a planet (Hatzes et al. 2003). 
\^ ' Aims. We aim to obtain a first direct detection of the stellar companion, to determine its current orbital position (for comparison with the 
spectroscopic and astrometric data), its infrared magnitude and, hence, mass. 

Methods. We use the Adaptive Optics camera CIAO at the Japanese 8m telescope Subaru on Mauna Kea, Hawaii, with the semi-transparent 
coronograph to block most of the light from the bright primary y Cep A, and to detect at the same time the faint companion B. In addition, we 
also used the IR camera CI Cass at the Calar Alto 3.5m telescope, Spain, to image y Cep A and B by adding up many very short integrations 



Oh' 
I 



O ' (without AO). 



Results, y Cep B is clearly detected on our CIAO and Q. Cass images. We use a photometric standard star to determine the magnitude of B 



C/2 after PSF subtraction in the Subaru image, and the magnitude difference between A and B in the Calar Alto images, and find an average value 
■ of A" = 7.3 + 0.2 mag. The separations and position angles between A and B are measured on 15 July 2006 and II and 12 Sept 2006, B is 
^ : slightly south of west of A. 

r Conclusions. By combining the radial velocity, astrometric, and imaging data, we have refined the binary orbit and determined the dynamical 
r^i masses of the two stars in the y Cep system, namely 1.40 ± 0.12 Mq for the primary and 0.409 ± 0.018 Mg for the secondary (consistent with 
^ being a M4 dwarf). We also determine the minimum mass of the sub-stellar companion to be M,, sin i = 1.60 ±0.13 Mj^p. 

Key words, instrumentation: adapive optics - binaries: spectroscopic - binaries: visual - planetary systems - star: individual: y Cep 



1. Introduction: y Cep 

The bright star y Cep (also HD 222404 or HIP 116727) lo- 
cated near the north celestial pole {a = 23''39'"20.8^ and 
(5 = -i-77°37'56.2" for J2000.0) is known as a single-lined spec- 
troscopic triple: The primary star has spectral type Kl 111-1 V, 
is visible to the naked eye {V = 3.2 mag) and is located at 
a distance of 13.8 pc. A low-mass stellar companion with an 
orbital period of several decades was discovered spectroscopi- 
cally by Campbell et al. (1988), but until now it has never been 
imaged directly (e.g. Hatzes et al. 2003, Mugrauer et al. 2006). 
Campbell et al. (1988) and Walker et al. (1992) discussed the 
evidence for additional radial velocity variations indicative of 
a very low-mass, possibly planetary companion. Hatzes et al. 
(2003) confirmed that there is indeed a third object in the sys- 
tem with a minimum mass of Mp sin i - 1 .7 Mjup, i.e., probably 
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a planet, with a semi-major axis of 2 AU and a 900-day eccen- 
tric orbit (e = 0.12). 

Several investigators have attempted to constrain the orbit 
as well as the mass of y Cep B (the stellar companion) on 
the basis of available spectroscopic data. The estimated pe- 
riod of the binary from these studies has ranged from about 
30 to 66 yr (Walker et al. 1992, Griffin et al. 2002). Most re- 
cently, Torres (2006) put together not only the spectroscopic 
data (dating back to 1902), but also Hipparcos and ground- 
based astrometric data (dating back to 1898) and obtained new 
results for the primary star (y Cep Aa), the secondary (B), and 
the sub-stellar companion (Ab): B orbits A in an eccentric or- 
bit (e = 0.4085 + 0.0065) with 66.8 + 1.4 year period and a 
semi-major axis of 19.02 + 0.64 AU; the primary has a mass 
of 1.18 + 0.11 M0, temperature of 4800 + 100 K, and an age 
of 6.6 Gyrs; the unseen stellar companion should have a mass 
of 0.362 + 0.022 Mq, i.e. spectral type about M4; the mini- 
mum mass of the lowest-mass component in the system is then 
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Fig. 1. Subaru CIAO K-band image of y Cep A and B, where 
A is below the semi-transparent coronograph. The total inte- 
gration time was 30 sec, obtained on 15 July 2006. 



Mp ■ sin; = 1.43 + 0.13 Mjup; the astrometric data yield an 
upper mass Umit of 17 Mjup (at the 3 cr confidence level), i.e. 
certainly sub-stellar, probably planetary; all values above are 
from Torres (2006). The predictions for 7 Cep B can be tested 
by direct imaging. 

In this paper, we present the first direct imaging detections 
of y Cep B (sect. 2), compare the parameters observed with the 
predictions (sect. 3), and newly determine the masses of y Cep 
Aa, Ab, and B (sect. 4). 

2. Imaging observations 

We observed y Cep first on 15 July 2006 with the 
Coronagraphic Imager for Adaptive Optics (CIAO) on the 
Japanese 8m telescope Subaru on Mauna Kea, Hawaii, in the 
K-band. The night was photometric with ~ 0.7 arc sec natu- 
ral seeing. We took ten images (always 30 times 0.1 sec each, 
i.e. 30 sec total integration time) through a semi-transparent 
coronographic mask (0.6 arc sec diameter, a chrome spot on a 
sapphire substrate) in order to block most of the Ught of y Cep 
A, but still to be able to detect B. 

The data were reduced in the normal way: Dark subtraction, 
flat-field devision, and adding up the individual images. The 
resulting image is shown in Fig. 1 . 

The semi-transparent coronograph has a transmission of 
between about 1% and 6%, thus still allowing the primary star 
to shine through. We subtracted the point spread function (PSF) 
of remaining light of y Cep A, see Fig. 2 for the result. 7 Cep 
B is clearly detected in Fig. 1 and 2. 

Because the throughput of the coronograph is not well 
known and because the star A may still be in the non-linear 
regime, we cannot derive the magnitude of B with respect to A 




Fig. 2. Subaru CIAO K-band image of 7 Cep A and B, now 
after subtraction of the PSF of y Cep A. The stellar companion 
7 Cep B is now clearly detected 0.87 arc sec west of A (slightly 
south of west). 



from the image itself. A photometric standard star is needed. 
We therefore observed the standard star FS 150 (9.94 + 0.01 
mag, UKIRT web page) in the same night (8x3 sec) to obtain 
the magnitude of 7 Cep B by aperture photometry (after PSF 
subtraction). The result for B is K = 7.32 ±0.11 mag. Hence, 
the magnitude difference between 7 Cep A and B is 6.28 ± 0.24 
mag (A has Ks = 1.04 ± 0.21 mag, 2MASS Catalog, Skrutskie 
et al. 2006). 

By Gaussian centering, we can determine the position of B 
relative to A, i.e. to measure the separation and position angle. 
For a companion (or candidate) to a very bright primary, it is 
very useful to observe with a semi-transparent coronograph, so 
that the primary star is still visible together with the companion 
(for separation measurement). For the astrometric analysis, we 
use the most recent pixel scale for CIAO, namely 21.3 + 0.1 
mas/pixel (Mayama et al. 2006) and the orientation of the de- 
tector being 0.8 1 ± 0.25 deg from north as measured by us from 
images of three sources in the V1686 Cyg region taken in the 
same night. 

The resulting separation and position angle are listed in 
Table 1. 

After our detection with Subaru we observed 7 Cep with the 
O Cass IR imager on the Calar Alto 3.5m telescope in Spain. 
Due to the fact that our instrument of choice, the AO system 
ALFA, was not functioning at the time, we observed 7 Cep 
in a backup programme by taking several thousand short (0. 1 
sec) integration images (speckle), which we then combined by 
shift-i-add, after background correction and flatfielding. We de- 
tected 7 Cep B on both 1 1 Sept (7391 images) and on 12 Sept 
2006 (18140) using the K' (K-prime) filter together with Br 7. 
In addition to measuring the separation between B and A from 
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Table 1. Astrometry and pliotometry for y Cep A and B 



Obs date 


Telescope 


PA 


O-C 


separation 


O-C 


^K 


HJD 


instrument 


n 


n 


[arc sec] 


[mas] 


[mag] 


2453932.046 


Subaru 8.3m CIAO 


256.91 + 0.27 


-0.20 


0.870 ± 0.005 


-0.005 


6.28 ± 0.24 


2453989.520 


CA 3.5m Q. Cass 


256.16 + 0.35 


+0.09 


0.891 +0.006 


+0.005 


6.28 + 0.13 


2453990.514 


CA 3.5m Q. Cass 


256.39 ± 0.34 


+0.33 


0.887 ± 0.005 


+0.001 


6.12 + 0.11 



Remark: Observed minus computed (O-C) residuals are explained in sect. 3. 
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Fig. 4. Orbit and orbital motion of y Cep B around A as already 
observed by our imaging within two month. The inlay shows 
our three imaging data points of B in orbit around A. 



the resulting images in each night, we determined also the mag- 
nitude difference relative to the primary, which was neither sat- 
urated nor in the non-linear regime of the detector. See Fig. 3 
for the final images of the two nights and Table 1 for the re- 
sults of our measurements. We used the Hipparcos multiples 
HIP 109474, HIP 91115, HIP 92961, HIP 104536, observed 
in the same two nights, for astrometric calibration of the pixel 
scale (38.681 +0.059mas/pixel) and the detector position angle 
(21.362 + 0.037°). 

Our measurements may indicate a small change in the po- 
sition angle of 0.63 + 0.36 deg between July and September 
2006, and a more significant increase in separation of 19 + 6 
mas, consistent with the predictions for orbital motion, see Fig. 
4. 

For the secondary star j Cep B, imaged for the first time 
here, we obtain a magnitude of = 7.3 + 0.2 mag, the average 
of the Subaru and the two Calar Alto data points. 

3. Comparison with predictions 

Our new measurements of the relative position of y Cep B with 
respect to the primary star allow us to test the predictions from 
the orbital solution for the binary produced recently by Torres 
(2006). That study combined all available measurements of the 



radial velocity of the primary from a variety of sources (includ- 
ing the high-precision studies that uncovered the sub-stellar 
companion, along with classical measurements from the litera- 
ture and from new observations) and used also the intermediate 
astrometry from the Hipparcos mission ('abscissa residuals') 
as well as ground-based astrometry reaching back more than 
a century. This enabled the scale and orientation of the binary 
orbit to be established for the first time, in a combined solu- 
tion that included also the perturbing effects of the sub-stellar 
companion. 

A comparison of our measurements with the published or- 
bit (Torres 2006) indicates systematic (observed minus com- 
puted) O-C residuals of ~ 8° in position angle and 50-60 
mas in separation (positive in both cases). These are many 
times larger than the measurement errors, although the pre- 
dicted positions themselves have uncertainties of about 3.1 de- 
grees in position angle and 34 mas in separation for the dates 
of our observations. These reflect the uncertainty in the or- 
bital elements as well as the estimated eiTor of the adopted pri- 
mary mass, which was external to the orbital solution of Torres 
(2006) and was established from stellar evolution considera- 
tions. Therefore, the residuals coiTespond to 1.5-2.5cr, which, 
while not overly large, may nevertheless indicate the need for 
some adjustment in the orbit given that they are systematic in 
nature. The new O-C values given in Table 1 (this paper) are 
computed for the new orbit (Table 2) after taking into account 
our three imaging observations. They are much smaller than 
in Torres (2006), i.e. a clear improvement as expected, which 
was the motivation for the imaging observations. The relative 
brightness of the secondary was predicted by Torres (2006) to 
be A/T ~ 6.4 mag, which is not far from the values we measure 
(see Table 1). 

4. Improved masses 

The assumption of a value for the primary mass in the Torres 
(2006) solution was necessary because of the lack of direct 
measurements of the relative position between j Cep A and B, 
since the secondary had never been imaged. With our present 
results this assumption can now be dropped, and this allows 
for a direct measurement of the dynamical masses of both stars 
with no need for external constraints. In order to update the or- 
bital solution and solve for the masses we have added to the 
data considered by ToiTes (2006) our 3 measurements of the 
separation and position angle of y Cep with their coiTespond- 
ing uncertainties. We have modified the procedures to include 
the relative semimajor axis of the binary, a^g, as an adjustable 
parameter in addition to the other variables considered in the 
original solution. The results are shown in Table 2, which in- 
cludes the most relevant elements of the inner and outer or- 
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Fig. 3. Cass IR speckle images of 7 Cep A and B, obtained on 1 1 Sept 2006 (left) and 12 Sept 2006 (right) with the Calar Alto 
3.5m telescope with ~ 12 to 30 min total integration time (through K' and Br y filters). 



bits and several derived quantities. The notation follows that of 
Torres (2006). Most of the orbital elements change very little 
(well within the errors). The most significant change is in the 
position angle of the ascending node (Oab). which increases by 
about 5°. This accounts for a large fraction of the residuals in 
P.A. noted above. The other significant change is in the scale of 
the orbit, which increases by about 6%. This in turn accounts 
for the residuals in the angular separation. The residuals from 
this new orbit are listed in Table 2, and are now all within lo", 
indicating a satisfactory fit. 

The resulting masses of the two stars are M^z. ' 1.40 ± 
0.12 Mo and Mb = 0.409 + 0.018 Mq. Both the primary 
and secondary are somewhat more massive than assumed by 
Torres (2006). Consequently the minimum mass of the sub- 
stellar companion also increases to Mp sin i = 1 .60 + 0.13 Mjup, 
a value intermediate between those proposed by Torres (2006) 
and Hatzes et al. (2003), and the upper mass limit from Torres 
(2006) changes from 17 Mjup to 19 Mj^p. The brightness of the 
stellar companion B is fully consistent with a single M4-type 
dwarf, neither multiple itself nor a white dwarf. 

Given the new parameters in Table 2, the critical semi- 
major axis for long-term stabihty of planets according to the 
Hohnan & Wiegert (1999) criterion is at only 3.85 ± 0.38 AU. 
The size of a planet-forming circumstellar disk in the close bi- 
nary is not larger than 3.4 AU according to the Pichardo et al. 
(2005) criteria. The planet candidate 7 Cep Ab orbits around 
Aa between 1.8 and 2.3 AU, which is inside the stable region, 
so that it is not neccessary to invoke migration after its forma- 
tion. 

Continued measurements of the relative position between 
7 Cep A and B over the coming years, which we intend to carry 
out, should further improve the stellar mass determinations and 
help estabUsh also the true mass of the sub-stellar companion. 



Table 2. Global orbital solution for y Cep 

Parameter and value 
Adjusted quantities from outer orbit (A+B) 

Pab (yr) 67.5 ± 1.4 

7 (km/s) -42.943 + 0.046 

ZA(km/s) 1.932 ±0.014 

Bab 0.4112 ± 0.0063 

WA(deg) 161.01 ±0.40 

rAB(yr) 1991.605 ±0.031 

a;(g(arcsec) 1.467 ±0.046 

a'l (mas) 332.4 ± 7.7 

lAB(deg) 119.3 ±1.0 

nAB(deg) 18.04 ±0.98 

Adjusted quantities from inner orbit (Aa+Ab) 

Pa (days) 902.9 + 3.5 

K,yjm/fi) 27.0 ±1.5 

gA 0.115 ±0.058 

aiAa (deg) 63 ± 27 

Ta (HJD-2,400,000) 53 147 ± 7 1 

Derived quantities 

ju;(masyr-') -63.86 ± 0.98 

Us (mas yr-') +150.76 ± 0.43 

n (mas) 72.69 ± 0.41 

Oab (AU) 20.18 ± 0.66 

MAa (M©) 1.40 ± 0.12 

Mb (Mq) 0.409 ± 0.018 

/(Mp)(10-''Mo)... 1.81 ±0.31 

MpSin;(Mjup) 1.60 ±0.13 

aAa-p(AU) 2.044 ±0.057 

Remark: The symbols have the same meaning as in Torres (2006). 
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